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Abstract 

We study circle valued maps and consider the persistence of the homology of their fibers. The 
outcome is a finite collection of computable invariants which answer the basic questions on persistence 
and in addition encode the topology of the source space and its relevant subspaces. Unlike persistence of 
real valued maps, circle valued maps enjoy a different class of invariants called Jordan cells in addition 
to bar codes. We establish a relation between the homology of the source space and of its relevant 
subspaces with these invariants and provide a new algorithm to compute these invariants from an input 
matrix that encodes a circle valued map on an input simplicial complex. 



'Department of Mathematics, The Ohio State University, Columbus, OH 43210,USA. Email: 
burgheleSmath . ohio-state . edu 

^Department of Computer Science and Engineering, The Ohio State University, Columbus, OH 43210, USA. Email: 

tamaldeyScse . ohio-state . edu 



1 Introduction 



Data analysis provides plenty of scenarios where one ends up with a nice space, most often a simplicial 
complex, a smooth manifold, or a stratified space equipped with a real valued or a circle valued map. The 
persistence theory, introduced in lTT3l . provides a great tool for analyzing real valued maps with the help 
of homology. Similar theory for circle valued maps has not yet been developed in the literature. The work 
in ll20l brings the concept of circle valued maps in the context of persistence by deriving a circle valued map 
for a given data using the existing persistence theory. In contrast, we develop a persistence theory for circle 
valued maps. 

One place where circle valued maps appear naturally is the area of dynamics of vector fields. Many 
dynamics are described by vector fields which admit a minimizing action (in mathematical terms a Lyapunov 
closed one form). Such actions can be interpreted as 1- cocycles which are intimately connected to circle 
valued maps as shown in (l). Consequently, a notion of persistence for circle valued maps also provides 
a notion of persistence for 1-cocycles which appear in some data analysis problems ETl l22l . In summary, 
persistence theory for circle valued maps promises to play the role for some vector fields as does the standard 
persistence theory for the scalar fields fl5] |6] IT31II91 . 

One of the main concepts of the persistence theoiy is the notion of bar codes irT9l -invaiiants that char- 
acterize a real valued map at the homology level. The angle (circle) valued maps, when characterized at 
homology level, require a new invariant called Jordan cells in addition to the refinement of the bar codes 
into four types. 

The standard persistence |[T3l [191 which we refer as sublevel persistence deals with the change in the 
homology of the sublevel sets which can not make sense for a circle valued map. However, the change in the 
homology of the level sets can be considered for both real and circle valued maps. The notion of persistence, 
when considered for the level sets of a real valued map (9] is referred here as level persistence. It refines the 
sublevel persistence. The zigzag persistence introduced in [4] provides complete invariants (bar codes) for 
level persistence of (tame) real valued maps. They are defined using representation theory for linear quivers. 

The change in homology of the level sets of a (tame) circle valued map is more complicated because 
of the return of the level to itself when one goes along the circle. It turns out that representation theory 
of cyclic quivers provides the complete invariants for persistence in the homology of the level sets of the 
circle valued maps. This notion of persistence is called here the persistence for circle valued maps and its 
invariants, bar codes and Jordan cells are shown to be effectively computable. 

Our results include a derivation of the homology for the source space and its relevant subspaces in terms 
of the invariants (Theorem 13 . 1 1 and [3 .2\ . The result also applies to real valued maps as they are special cases 
of the circle valued maps. This leads to a result (Corollary 13.41 ) which to our knowledge has not yet appeared 
in the literature Q. A number of other topological results which can not be derived from any of the previously 
defined persistence theories are described in [3 j providing additional motivation for this work. 

After developing the results on invariants, we propose a new algorithm to compute the bar codes and 
Jordan cells. For a simplicial complex, the entire computation can be done by manipulating the original 
matrix that encodes the input complex and the map. The algorithm first builds a block matrix from the orig- 
inal incidence matrix which encodes linear maps induced in homology among regular and critical level sets, 
more precisely the quiver representations p r described in section @] Next, it iteratively reduces this new ma- 
trix eliminating and hence computing the bar codes. The resulting matrix which is invertible can be further 
processed to Jordan canonical form [101, providing Jordan cells. The algorithm for zigzag persistence [@] 
when applied to what we refer in section 3 as the infinite cyclic covering map f can compute bar codes 
but not Jordan cells. In contrast, our method can compute the bar codes and Jordan cells simultaneously by 

it was brought to our attention by David Cohen-Steiner that the extended persistence proposed in (6) allows similar connections 
between homology of source spaces and persistence. 
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manipulating matrices and can also be used as an alternative to compute the bar codes in zig-zag persistence. 

Notations. We list here some of the notations that are used throughout. 

• For rth homology group of a topological space X under an a priori fixed field k, we write H r (X) 
instead of H r (X; k). 

• For a map / : X — > Y and K CY we write Xk ■= 

• We use Z>o and Z>o for non-negative and positive integers respectively. 

• In our exposition, we need to use open, semi-open, and closed intervals denoted as (a, b), (a, 6] or 
[a, b), and [a, b] respectively. To denote an interval, in general, we use the notation {a, b} where "{" 
stands for either "[" or "(". 

• For a linear map a : V —> W between two vector spaces we write : 

kera := {v 6 V | a(v) = 0}, img a := {w € a(V) C W}, coker a := W/a(V). 

• A matrix A is said to be in column echelon form if all zero columns, if any, are on the right to nonzero 
ones and the leading entry (the first nonzero number from below) of a nonzero column is always 
strictly below of the leading entry of the next column. Similarly, A is said to be in row echelon form 
if all zero rows, if any, are below nonzero ones and the leading entry (the first nonzero number from 
the right) of a nonzero row is always strictly to the right of the leading entry of the row below it. 

If A is an m x n matrix (m rows and n columns), there exist an invertible n x n matrix R(A) and 
an invertible m x m matrix L(A) so that A ■ R(A) is in column echelon form and L(A) ■ A is in row 
echelon form. Algorithms for deriving the column and row echelon form can be found in standard 
books on linear algebra. 

2 Definitions and background 

We begin with the technical definition of tameness of a map. 

For a continuous map / : X — > Y between two topological spaces X and Y, let Xjj = f~ l (U) for 
U C Y. When U = y is a single point, the set X y is called a. fiber over y and is also commonly known as the 
level set of y. We call the continuous map / : X — > Y good if every y 6 7 has a contractible neighborhood 
U so that the inclusion X y — > Xjj is a homotopy equivalence. The continuous map / : X — > Y is afibration 
if each y E Y has a neighborhood U so that the maps / : Xjj — >• U and pr : X y x U — > U are fiber wise 
homotopy equivalent. This means that there exist continuous maps / : Xjj — > X y x U with pr\u ■ l\u = f\u 
which, when restricted to the fiber for any z G U, are homotopy equivalences. In particular, / is good. 

Definition 2.1 A proper continuous map f : X — > Y is tame if it is good, and for some discrete closed 
subset S C Y, the restriction f : X \ / _1 (S I ) — > Y \ S is afibration. The points in S CY which prevent f 
to be afibration are called critical values. 

If Y = R and X is compact or Y = S 1 , Hthen the set of critical values is finite, say s± < S2 < ■ ■ ■ Sfc. 
The fibers above them, X Si , are referred to as singular fibers . All other fibers are called regular. In the case 
of S 1 , Si can be taken as angles and we can assume that < Sj < 2tt. Clearly, for the open interval Sj) 
the map / : / -1 (sj_i, Si) — >• (sj-i, Sj) is a fibration which implies that all fibers over angles in Sj) 
are homotopy equivalent with a fixed regular fiber, say X tl , with t,- L G (sj_i, si). 

2 since the map / is proper and S 1 compact, so is X 
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X ti a, 6j In particular, there exist maps cjj : — > X s . and &j : — > 

X Si , unique up to homotopy defined as follows: If t{ and U + i are 
contained in Ui C Y where the inclusion X Si C Xu i is a homotopy 
equivalence with a homotopy inverse n : Xu i — > X Si , then a, and b{ 
are the restrictions of to X ti and respectively. If not, in view 
of the tameness of /, one can find t\ and t' i+1 in U so that X ti and X ti+1 are homotopy equivalent to X t / 
and X t / i respectively and compose the restrictions of with these homotopy equivalences. These maps 
determine homotopically / : X — > Y, when Y = M. or S 1 . For simplicity in writing, when F = Kwe put 
tfc+i € (sfe, oo) and ii € (—00, Si) and when Y = S 1 we put t^ + \ = t\ £ (s&, si + 27r). All scalar or circle 
valued simplicial maps on a simplicial complex, and all smooth maps with generic isolated critical points 
on a smooth manifold or stratified space are tame. In particular, Morse maps are tame. 

2.1 Persistence and invariants for real valued maps 

Since our goal is to extend the notion of persistence from real valued maps to circle valued maps, we first 
summarize the questions that the persistence answers when applied to real valued maps, and then develop a 
notion of persistence for circle valued maps which can answer similar questions and more. We fix a field k 
and write H r (X) to denote the homology vector space of X in dimension r with coefficients in a field k. 

Sublevel persistence. The persistent homology introduced in |[T3l and further developed in |fl9l is con- 
cerned with the following questions: 

Ql. Does the class x G H r (X^_ oo t ^) originate in H r (X^_ oot i^) for t" < tl Does the class x € 
#r(^(-oo,t]) vanish in H^X^^) for t < t'l 

Q2. What are the smallest t' and largest t" such that this happens? 

This information is contained in the inclusion induced linear maps H r {X^_ (X> ^) — > i?r(^(-oo,t']) where 
t' > t and is known as persistence. Since the involved subspaces are sublevel sets, we refer to this persis- 
tence as sublevel persistence. When / is tame, the persistence for each r = 0, 1, • • • dim X, is determined by 
a finite collection of invariants referred to as bar codes |[T9l . For sublevel persistence the bar codes are a col- 
lection of closed intervals of the form [s, s'] or [s, 00) with s, s' being the critical values of /. From these bar 
codes one can derive the Betti numbers of ^(_oo,a]) tne dimension of img(H r (X(_ 00 t j) — > H^X/jo^)) 
and get the answers to questions Ql and Q2. For example, the number of r-bar codes which contain the 
interval [a, b] is the dimension of img(H r (X(_ 00 a j) — > H r (X(_ 00 b ])). The number of r-bar codes which 
identify to the interval [a, b] is the maximal number of linearly independent homology classes bom exactly 
in -X"(_ 0Oi(l ] but not before and die exactly in H r (X_ 00 ^) but not before. 

Level persistence. Instead of sublevels, if we use levels, we obtain what we call level persistence. The 
level persistence was first considered in |Q but was better understood computationally when the zigzag 
persistence was introduced in [@]. Level persistence is concerned with the homology of the fibers H r (X t ) 
and addresses questions of the following type. 

Ql. Does the image of x € H r (X t ) vanish in H r (Xu t n), where t' > t or in H r (Xun ^1), where t" < tl 

Q2. Can x be detected in H r (X t >) where t' > t or in H r (X t n) where t" < tl The precise meaning of 
detection is explained below. 

Q3. What are the smallest t' and the largest t" for the answers to Ql and Q2 to be affirmative? 
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To answer such questions one needs information about the following inclusion induced linear maps: 

H r (X t ) ^ H r (X m ) ^ H r (X t ,). 

The level persistence is the information provided by this collection of vector spaces and linear maps for all 
t,t'. 

We say that x € H r {X t ) is dead in iJ r (Xr tjt /i), t' > t, if its image by H r (X t ) — > H r (X\ t t n) vanishes. 
Similarly, x is dead in H r {X^h ^ ) , t" < t, if its image by H r (X t ) — >• H r (X^ t n ^) vanishes. 

We say that x € H r (X t ) is detected in H r (X t /), t' > t, (resp. t" < t), if its image in H r {Xu^n ) 
(resp. in H r (Xur/£) is nonzero and is contained in the image of H r {X t i) — > H r (X^ t ^) (resp. H r (X t ») — > 
H r (X[ t « j])). In Figure [Q the class consisting of the sum of two circles at level t is not detected on the 
right, but is detected at all levels on the left up to (but not including) the level t'. In case of a tame map the 
collection of the vector spaces and linear maps is determined up to coherent isomorphisms by a collection of 
invariants called bar codes for level persistence which are intervals of the form [s, s'], (s, s'), (s, s'], [s, s') 
with s,s' critical values as opposed to the bar codes for sublevel persistence which are intervals of the 
form [s, s'], [s, oo) with s, s' critical values. These bar codes are called invariants because two tame maps 
/ : X — > M and g : Y — > M which are fiber wise homotopy equivalent have the same associated bar codes. 
In the case of level persistence the open end of an interval signifies the death of a homology class at that 
end (left or right) whereas a closed end signifies that a homology class cannot be detected beyond this level 
(left or right). In the case of the sublevel persistence the left end signifies birth while the right death. Level 
persistence provides considerably more information than the sub level persistence. The bar codes of the sub 
level persistence can be recovered from the ones of level persistence. Precisely a level bar code [s, s'} gives a 
sublevel bar code [s, oo) and a level bar code [s, s') gives a sublevel bar code [s, s']; the sublevel persistence 
does not see any of the level bar codes (s, s') or (s, a']. It turns out that the bar codes of the level persistence 
can also be recovered from the bar codes of the sub level persistence of / and additional maps canonically 
associated to /. 

In Figure [T] we indicate the bar codes both for sub level and level persistence H for some simple map 
/ : X — )• R in order to illustrate their differences. The space X is a tube open on one end and / is the height 
function laid horizontally. 

D 



level persistence H ^ 



sub-level persistence |_ 

f t 
Figure 1: Bar codes for level and sub-level persistence. 



3 the white circles indicate open ends and the dark circles indicate closed ends 
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3 Persistence for circle valued maps 



Let / : X — > S 1 be a circle valued map. The sublevel persistence for such a map cannot be defined since 
circularity in values prevents defining sub-levels. Even level persistence cannot be defined as per se since 
the intervals may repeat over values. To overcome this difficulty we associate the infinite cyclic covering 
map / : X — > R for /. It is defined by the commutative diagram: 



X 



X 



f 



f 



The map p : R — > S 1 is the universal covering of the circle (the map which assigns to the number t G K the 
angle 6 = t(mod 2ir) and tp is the pull back of p by the map / which is an infinite cyclic covering. Notice 
that if p(t) = 9 then X t and Xq are identified by ip. If x G H r (Xg) = H r (X t ), p(t) = 9, the questions Ql, 
Q2, Q3 for / and X can be formulated in terms of the level persistence for / and X. 

Suppose that x G H r (Xt) = H r (Xg) is detected in H r (X t >) for some t' > t + 2-k. Then, in some 
sense, x returns to H r (Xe) going along the circle 8 1 one or more times. When this happens, the class x 
may change in some respect . This gives rise to new questions that were not encountered in sublevel or level 
persistence. 

Q4. When x G H r (Xg) returns, how does the "returned class" compare with the original class xl It may 
disappear after going along the circle a number of times, or it might never disappear and if so how 
does this class change after its return. 

To answer Q1-Q4 one has to record information about H r (Xg) — > H r {X\Q^iA <— H r (Xgi) for any pair of 
angles 9 and 9' which differ by at most 2ir. This information is referred to as the persistence for the circle 
valued map f. 

When / is tame, this is again completely determined up to coherent isomorphisms by a finite collection 
of invariants. However, unlike sublevel and level persistence for real valued maps, the invariants include 
structures other than bar codes called Jordan cells. Specifically, for any r = 0, 1, • • • , dim(X) we have two 
types of invariants: 

• bar codes: intervals with ends s,s' < s < 2-k, s < s' < oo, that are closed or open at s or s', 
precisely of one of the forms [s, s'], (s, s'], [s, s'), and (s, s'). These intervals can be geometerized as 
"spirals" with equations in CD- For any interval {s, s'} the spiral is the plane curve (see Figure |3]in 
section |4]) 



x{6) = (9 + 1- s)cos9 
y (9) = (# + l-s)sin# 



with 9 G {s, s'}. 



(1) 



Jordan cells. A Jordan cell is a pair (A, k), A G k\0, k G Z>o, where k denotes the algebraic closure 
of the field k. It corresponds to a k x k matrix of the form 



/A 1 0... 0\ 
A 1... 







A 




f 

V 



(2) 
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• r-invariants. Given a tame map / : X — > § , the collection of bar codes and Jordan cells for each 
dimension r G {0, 1, 2, • • • dimX} constitute the r-invariants of the map /. 

We will define all of the above items in the next section using quiver representations. 

The bar codes for / can be inferred from / : X\ a n — > R with [a, b] being any large enough interval. 
Specifically, the bar codes of / : X — >• S 1 are among the ones of / : Xy a ^ — > R for (b — a) being at most 
sup 9 dim H t (Xq). 

The Jordan cells can not be derived from / : X — > R or any of its truncations / : X\ a u — > R unless 
additional information, like the deck transformation of X, is provided. The end points of any bar code for / 
correspond to critical angles, that is, s and s' (mod 27r) of a bar code interval {s, s'} are critical angles for 
/. One can recover the following information from the bar codes and Jordan cells: 

1. The Betti numbers of each fiber, 

2. The Betti numbers of the source space X, and 

3. The dimension of the kernel and the image of the linear map induced in homology by the inclusion 
Xq C X as well as other additional topological invariants not discussed here 0. 

Theorems 13.11 and 13.21 make the above statement precise. Let B be a bar code described by a spiral 
in (Q]) and 9 be any angle. Let nQ(B) denote the cardinality of the intersection of the spiral with the ray 
originating at the origin and making an angle 9 with the x-axis. For the Jordan cell J = (A, k), let n( J) = k 
and A (J) = A. Furthermore, let B r and J r denote the set of bar codes and Jordan cells for r-dimensional 
homology. We have the following results. 

Theorem 3.1 dimH r (X e ) = EeeBr n o( B ) + Ejgj r n ( J )- 

Theorem 3.2 dimH r (X) = #{B G fi r |both ends closed} + #{B G # r _i|both ends open} + #{J € 
J r \X(J) = 1} + #{J G J r -i\\{J) = 1}. 

Using the same arguments as in the proof of the above Theorems one can derive: 

Proposition 3.3 dimimg(H r (X e ) H r (X)) = #{B G fi r |n e (B) / and both ends closed} + #{J G 
Jr\\ = 1} 

A real valued tame map / : X — > R can be regarded as a circle valued tame map /' : X — >■ S 1 
by identifying R to (0, 2-k) with critical values t\, ■ ■ ■ ,t m becoming the critical angles 0\, ■ ■ ■ ,9 m where 
9{ = 2 arctan t j + 7r. The map /' in this case will not have any Jordan cells and the bar codes will be the 
same as level persistence bar codes. We have the following corollary: 

Corollary 3.4 dimH r (Xg) = J^BeBr n s( B ) an ^ 

dimi^AT) = #{B G i3 r |both ends closed} + #{B G # r _i|both ends open}. 

Theorem 13.11 is quite intuitive and is in analogy with the derived results for sublevel and level persis- 
tence El QjO. Theorem l3.2l is more subtle. Its counterpart for real valued function (Corollary 13.41 has not yet 
appeared in the literature though a related result for homology of source space can be derived from extended 
persistence O. The proofs of these results require the definition of the bar codes and Jordan cells which 
appear in the next section. The proofs are sketched in section [5] 

The Questions Q1-Q3 can be answered using the bar codes. The question Q4 about returned homology 
can be answered using the bar codes and Jordan cells. 
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r-invariants 


circle 1: 1 times around 1,-3 times around 2, -2 times around 3 
circle 2: 1 times around 1, 4 times around 2, 1 time around 3 
circle 3: 2 times around 1, 2 times around 2, 2 times around 3 


dimension 


bar codes 


Jordan cells 







(1,1) 


1 


(0 6 ,0i + 2tt] 

[02,03] 
(04,0b) 


(3,2) 



Figure 2: Example of r-invariants for a circle valued map 



Figure |2] indicates a tame map / : X — > S 1 and the corresponding invariants, bar codes, and Jordan 
cells. The space X is obtained from Y in the figure by identifying its right end Y\ (a union of three circles) 
to the left end Yq (again a union of three circles ) following the map (j) : Y\ — > Yq. The map / : X — s> S 1 
is induced by the projection of Y on the interval [0, 2ir]. We have Hi(Y±) = H\{Yq) and 6 

induces a linear map in 1 -homology represented by the matrix 

n i 2 \ 

-3 4 2 . 
V-2 1 2/ 

The first generator (circle 1) of H\{X.2 V ) is dead in H\{X\a 2t ]) for < 9q but not for ^ G (^6, 27r] and is 
detected in Hi^X^^+e) for < < 9\ but not for 9 > 9\. It generates a bar code 27r + ^i]. The other 
two (circle 2 and 3) never die and provide a Jordan cell (3, 2). In Appendix we show how our algorithm can 
be used to compute the bar codes and Jordan cells for the above example. 



4 Representation theory and its invariants 

The invariants for the circle valued map are derived from the representation theory of quivers. The quivers 
are directed graphs. The representation theory of simple quivers such as paths with directed edges was 
described by Gabriel ifTTI and is at the heart of the derivation of the invariants for zigzag and then level 
persistence in |4). For circle valued maps, one needs representation theory for circle graphs with directed 
edges. This theory appears in the work of Nazarova lUTI . and Donovan and Ruth-Freislich lfl2l . 

Each circle is oriented counterclockwise and represents a 1-dimensional homology class; "k times (—k times) around the 
circle" means " going around k times counter clockwise (clockwise respectively)". 
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Let G 2m be a directed graph with 2m vertices, xi, x\, ■ ■ ■ x%. 



i 




x 2 
ax] 



Its underlying undirected graph is a simple cycle. The directed 
edges in G 2m are of two types: forward ai : x 2 i^\ — > x 2 i, 1 < i < 
m, and backward bi : x 2 i + \ — » x 2 i, 1 < i < m— 1, b m : x\ — > a^m- 



We think of this graph as being located on the unit circle cen- 
tered at the origin o in the plane. 



X 2m -2 



X 2m -1 



A representation p on G 2m is an assignment of a vector space 
V x to each vertex x and a linear map £ e : V x — >• for each oriented 
edge e = {x, y}. Two representations p and // are isomorphic if 
for each vertex x there exists an isomorphism from the vector space 
V x of p to the vector space V x of p' , and these isomorphisms com- 
mute with the linear maps V x — > V y and V x — > V' A non-trivial 



representation assigns at least one vector space which is not zero-dimensional. A representation is indecom- 
posable if it is not isomorphic to the sum of two nontrivial representations. 

Given two representations p and p', their sum p © p' is a representation whose vector spaces are the 
direct sums V x © V x related by linear maps that are the direct sums t e © £' e . It is not hard to observe that each 
representation has a decomposition as a sum of indecomposable representations unique up to isomorphisms. 

We provide a description of indecomposable representations of the quiver G 2m . For any triple of integers 
{i,j, k}, 1 < i, j < m, k > 0, one may have any of the four representations, p k), j]; k), 

j); k) , and p k) defined below. For any Jordan cell (A, A;) one has the representation p J (\, k) 

defined below. The exponents / and J indicate that these representations are associated with a bar code 
(interval) or a Jordan cell respectively and hence we call them bar code and Jordan cell representations. 

• Bar code representation /f/({i, j}; k): Suppose that the evenly indexed vertices {x 2 , x^, • • • x 2m } of 
G 2m which are the targets of the directed arrows correspond to the angles < s\ < s 2 < ■ ■ ■ < s m < 
2tt. Draw the spiral curve given by £Q for the interval {si, sj + 2kir}; refer to Figure [3] 

For each Xi, let {ej, e|, ■ ■ ■ } denote the ordered set (possibly empty) of intersection points of the ray 
oxi with the spiral. While considering these intersections, it is important to realize that the point 
(x(si),y(si)) (resp. (x(sj + 2kir),y(sj + 2kir))) does not belong to the spiral (0Q) if {i, j} is open at 
i (resp. j). For example, in Figure [3] the last circle on the ray ox 2 j is not on the spiral since [i, j) in 
p T ([i,j)\ 2) is open at right. 

Let V Xi denote the vector space generated by the base {ej,ef, ■ ■ ■ }. Furthermore, let Oj : V X2i _ 1 — > 
V X2i and Pi : V X2i+1 — > V X2i be the linear maps denned on bases and extended by linearity as follows: 
assign the vector S V Xi to e^^i if e 2i is an adjacent intersection point to the points e 2i ±i on the 
spiral. If e 2 \ does not exist, assign zero to e\ i±l . If e| i _ ) _ 1 do not go to zero, h has to be I, I — 1, or 
I + 1. The construction above provides a representation on G 2m which is indecomposable. Once the 
angles Sj are associated to the vertices x 2 i one can also think of these representations p 1 ({i, j}; k) as 
the bar codes [s^ Sj + 2kn], (sj, Sj + 2/c7r], [sj, Sj + 2A;7r), and (sj, sj + 2kn). 

• Jordan cell representation p J (A, k): Assign the vector space with the base {e±, e 2 , ■ ■ ■ , e^} to each Xi 
and take all linear maps a.; but one (say a.\) and the identity. The linear map ot\ is given by the 
Jordan matrix defined by (A, k) in (2). Again this representation is indecomposable. 

It follows from the work of lfT2l [T71 that when k is algebraically closed^, the bar code and Jordan cell 
representations are all and only indecomposable representations of the quiver G 2m . The collection of all bar 
code and Jordan cell representations of a representation p constitutes its invariants. 

5 when « is not algebraically closed Jordan cells have to be replaced by conjugate classes of indecomposable (not conjugated to 
a direct sum of matrices) matrices with entries in ft. 
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Figure 3: The spiral for [sj, sj + 47r). 



Now, consider the representation p on the graph G 2m given by the vector spaces V21-1 := V X2i _ 1 , 
V^ 2i and the linear maps a-i and To such a representation p, we associate a map M p : © 1<i<m V-^-i 
0l<i<m which is represented by a block matrix also denoted as M p : 



( oci 




-Pi o 

a 2 -fa 



\ 




a m -i -p m -i 

\ fim ■ ■ ■ ■ • • C^m j 

For this representation we define its dimension characteristic as the 2m-tuple of positive integers 

dim(p) = (n l5 r\ ■ ■ ■ n m , r m ) 

with rii = dimV X2i _ 1 and rj = dimV X2i and denote by ker(p) := kerM p and coker(/)) = cokerM p . For 
the sum of two such representations p = p\ © p 2 we have: 

Proposition 4.1 

7. dim(/?i © p 2 ) = dim(pi) + dim(/3 2 ), 

2. dimker(/?i © p 2 ) = dimker( / oi) + dimker^). 

3. dimcoker(pi © p 2 ) = dimcoker(/?i) + dimcoker^)- 

The description of a bar code representation permits explicit calculations. 

Proposition 4.2 

7. Ifi < j then 
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(a) dim p 1 k) is given by: 

ni = k + 1 if {i + 1) < I < j and k otherwise, 
ri = k + lifi < I < j and k otherwise 

(b) dim p 1 k) is given by: 

rii = k + 1 if {i + 1) < I < j and k otherwise, 
ri = k + 1 if (i + 1) < I < j and k otherwise, 

(c) dim p 1 ( [i , j ) ; k) is given by: 

ni = k + 1 if '{i + 1) < I < j and k otherwise, 
r\ = k + lifi < I < {j — 1) arcc? A; otherwise, 

( d) dim p 1 ((i,j);k) is given by: 

ni = k + 1 if {i + 1) < I < j and k otherwise, 

ri = k + 1 if (i + 1) < I < (j — 1) aTifif otherwise 

2. Ifi>j then similar statements hold. 

(a) dim // ( [i, j] ; k) is given by: 

n\ = k if (j + 1) < I < i and k + 1 otherwise; 

ri = k if (j + 1) < I < (i — l)j and k + 1 otherwise 

(b) dim p 7 ((i, j]; k) is given by: 

ni = k if (j + 1) < I < i and k + 1 otherwise. 
T[ = k if (j + 1) < I < i and k+1 otherwise, 

(c) dim p 1 ( [i , j ) ; k ) is given by: 

ni = k if (j + 1) <l<i and k+1 otherwise; 
ri = k if j < I < (i — 1) and k + 1 otherwise, 

(d) dim j); fc) z's given by: 

ni = k if [j + 1) < / < i and k+1 otherwise; 
ri = k if j < I < i and k + 1 otherwise. 

3. dim p J (X, k) is given by ni = ri = k 

Proposition 4.3 

1. dimker j]; k) = 0, dimcoker// ([i, j]; k) = 1, 

2. dimker p 1 ([i, j); k) = 0, dimcokeiy/Qi, j); k) = 0, 

3. dimker p J ((i,i]; k) = 0, dimcoker//((i, j]; k) = 0, 

4. dimker k) = 1, dimcoker// ((i, j); k) =0, 

5. dimker p J (\, k) = fresp. 1) if 1 (resp. 1), 

6. dimcoker/9 J (A, k) = (resp. 1) if 1 (resp. 1). 
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Observation 4.4 A representation p has no indecomposable components of type p in its decomposition iff 
all linear maps a[s and (3[s are isomorphisms. For such a representation, starting with an index i, consider 
the linear isomorphism 

Ti = /3" 1 • at ■ j3~\ ■ cti-i ■ ■ ■ p^ 1 ■ a 2 • /Jf 1 • «i • /3" 1 • a m ■ • a m -i • • • /3^ 1 • a i+ i. 

The Jordan canonical form H10\l of the isomorphism Ti is independent of i and is a block diagonal matrix 
with the diagonal consisting of Jordan cells (A, k)s. Clearly, p is the direct sum ofp J (X, k)s, the Jordan cell 
representations of p. 

Definition 4.5 (r-invariants.) Let f be a circle valued tame map defined on a topological space X. For f 
with m critical angles < s\ < S2, • • • s m < 2tt, consider the quiver G2m with the vertices xn identified 
with the angles s, and the vertices xi%—\ identified with the angles ti that satisfy < t\ < si < t% < 

For any r, consider the representation p r ofG^m with V Xi = H r (X Xi ) and the linear maps oiiS and fas 
induced in the r-homology by maps ai : X X2i _ 1 — > X X2i and bi : X X2i+l — > X X2i described in section\2\ 
The bar code and Jordan cell representations of p r are independent of the choice oftiS and are collectively 
referred as the r-invariants of f. 



5 Proof of the main results 

The Figure |2] and the bar codes listed below suggest why a semi-closed (one end open and the other closed) 
bar code does not contribute to the homology of the total space X and why a closed bar code (both ends 
closed) in B r contributes one unit while an open (both end open) bar code in £> r _i contributes one unit to the 
H r {X). Indeed, in our example, a semi-closed bar code in B\ adds to the total space a cone over S 1 , which 
is a contractible space. It gets glued to the total space along a generator of the cone (a segment connecting 
the apex to S 1 ), again a contractible space. A closed bar code in B\ adds a cylinder of S 1 whose H\ has 
dimension 1. It gets glued to the total space along a generator of the cylinder (a segment connecting the 
same point on the two copies of S 1 ), again a contractible space. An open bar code in B\ adds the suspension 
over S , topologically a 2-sphere which gets glued along a meridian, a contractible space. This contributes 
a dimension to H2. 

The lack of contribution of a Jordan cell with A 7^ 1 as well as the contribution of one unit of a Jordan 
cell in J r with A = 1 to both r and r + 1 dimensional homology of the total space should not be a surprise 
for the reader familiar with the calculation of the homology of mapping torus. 

Below we explain rigorously but schematically the arguments for the proof of Theorems 13.11 13.21 and 
Corollary EH 

The proof of Theorem 13.11 is a consequence of Propositions 14.11 and 14.21 The proof of Theorem 13.21 
proceeds along the following lines. 

First observe that, up to homotopy, the space X can be regarded as the iterated mapping torus T de- 
scribed below. Consider the collection of spaces and continuous maps: 

V V b <-) = bm D ai K V M D a2 K V V b m~l D O-m v 

-&m — Ao < rtl 1 Ai i K2 > A2 • • ■ A m _i i K m > A m 

with Ri := X ti and Xj := X Si and denote by T = T(ol\ • • • a m ; /3± ■ • ■ /3 m ) the space obtained from the 
disjoint union 

( y J?* x[0,l])U( □ Xi) 

l<i<m l<i<m 

by identifying Ri x {1} to Xi by a, and Ri x {0} to AV-i by /3j_i. Denote by f T : T — > [0, m] where 
f T : Ri x [0, !]—>[« — 1, i] is the projection on [0, 1] followed by the translation of [0, 1] to [i — 1, i\. This 
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map is a homotopical reconstruction of / : X —> S 1 provided that, with the choice of angles ij, and maps 
a, described in section|2j X, := f^ 1 (si), Ri := f~ l (ti). 
Let T 7 ' denote the space obtained from the disjoint union 

( y i2iX(e,l])U( □ X,) 

l<i<m l<j<m 

by identifying i?j x {1} to Xi by aj, and V" denote the space obtained from the disjoint union 

( □ ^x[0,l-e)U( □ X,) 



KKm 



Ki<m 



by identifying Ri x {0} to by 

Let K = Ul<i< m R i and X = Ul<i<m X i- Then > one has: 

1. r = v'uv", 

2. V' n V" = (Ui<i< m Ri x (e, 1 - e)) U ^, and 

3. the inclusions (Ui<i< m #i x { 1 / 2 }) u <*' C V'V\V" as well as the obvious inclusions # C P'and C 
V" are homotopy equivalences. 



The Mayer- Vietoris long exact sequence leads to the diagram 



H r (K) 



Mo 



H r {X) 



pri 



{Id-Id) 



H r+1 (T) © H r (X) £T r (^) © H r {Jfy-^ H r (T) 



H r {X) 



Id 



H r {X) 



Here A denotes the diagonal, in 2 the inclusion on the second component, pr\ the projection on the first 
component, V the linear map induced in homology by the inclusion X C T, and M Pr the map given by the 
matrix 



/ «i -ft o 
a r 2 



\ 




(3) 



Cl 

\-/3 r a r / 

\ ™ • • • • • • ■ ■ ■ ■ ■ ■ "m / 

with a[ : H r {Ri) — > H r (Xj) and : H r (Ri + \) — > H r (Xi) induced by the maps ctj and and iV defined 
by 



where a r and /3 r are the matrices 



a r Id 
-P r Id 



(a\ \ 

a r 2 

\0 ... a r m _J 



12 







p\ ... \ 
o pi ... 




M. 



Pr 



M, 



>H r {K) 



> iJ r (Af) -> H r {T) -> ff r _i(K) 



».ff r _i(#) -> ••• 



(4) 



from which we derive the short exact sequence 



-> coker(/9 r ) -> H r (T) ->■ ker(/3 r _i) ->■ 



(5) 



and then 



H r (T) = coker/3 r © ker p r _i 



(6) 



Theorem l3.2l follows from Propositions 14. 1 1 1431 and the equation © above. A specified decomposition 
of p r and p r _i into indecomposable representations and a splitting in the sequence © provide specified 
elements in H t {Xq) and H r (T) which can be compared. This leads to the verification of Proposition [33] 

6 Algorithm 

Given a circle valued tame map / : X — > § , we now present an algorithm to compute the bar codes and the 
Jordan cells when X is a finite simplicial complex, and / is generic and linear. This makes the map tame. 
Genericity means that / is injective on vertices. To explain linearity we recall that, for any simplex a £ X, 
the restriction f\ a admits liftings / : a — >• R, i.e. / is a continuous map which satisfies p- f = f\ a . The map 
/ : X — > S 1 is called linear if for any simplex a, at least one of the liftings (and then any other) is linear. 

Our algorithm takes the simplicial complex X equipped with the map / as input and, for any r, 
computes the matrix M Pr of the representation p r for /. This requires recognizing the critical values 
s\i S2, ■ ■ ■ s m 6 S 1 of /, and for conveniently chosen regular values t±, £2, ■ ■ ■ t m € S , determining the 
vector spaces V21-1 = H r (X ti ), V^i = H r (X Si ) with the linear maps a» and /3j as matrices. We consider 
the block matrix M Pr : 1<i<m Vn-\ — >• 0i<i< m Vu described in the previous section. 

We compute the bar codes from the block matrix M Pr first, and then the Jordan cells. The algorithm 
consists of three steps. We describe the first and second steps in sufficient details. The third step is a routine 
application of Observation 4. 1 and is accomplished by standard algorithms in linear algebra (reduction of 
the matrix to the canonical Jordan form). 

• Step 1. Compute the matrices 014, Pi that constitute the matrix M Pr of the representation p r . 

• Step 2. Process the matrix of M Pr to derive the bar codes ending up with a representation p' r whose 
all a-s and fys are invertible matrices. 

• Step 3 Compute the Jordan cells of p r from the representation p' r . 
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Step 1. In Step 1 we begin with the incidence matrix of the input simplicial complex X equipped with the 
map / : X — > S 1 . Let the angles < s\ < S2 ■ ■ ■ s m < 2ir be the critical values of /. Choose a collection 
of regular angles < t\ < t2 • • • t m < 2ir with ti < Sj < ti + i < Sj+i. Consider a canonical subdivision 
of X into a cell complex so that Xu^ t i, and X tt are subdivided into subcomplexes Ri and Xj as follows. 
For any open simplex a we associate the open cells : 

1. a(i) := a n with dim(cr(i)) = dime" — 1 if the intersection is nonempty 

2. <r(i) := a Pi -XV^t ) with dim cr(i) = dimcr if the intersection is nonempty. 

The cells of Xi are exactly of the form a(i) and their incidences are given as I(a(i),T(i)) = I{ct,t) 
where I(a, r) = 0, +1, or — 1 depending on whether r is a coface of a and whether their orientations 
match or not. The cells of Ri consist of cells of Xi, Xi + \, and all cells of the form cr(i). The incidences 
are given as I (a (i) , r (i)) = J(cr, r), I(a(i),a(i)) = 1, and + 1), = —1- All other incidences 
are zero. Assume that we are given a total order for the simplices of X that is compatible with / and 
also the incidence relations. This induces a total order for the cells in Xi and Xi + i and also the cells in 
R' t = Ri \ Xi U Xi + i for any 1 < i < m with X m+ \ := X±. Impose a total order on Ri by juxtaposing the 
total orders of Xi, -Xj+i, and R[ in this sequence. Clearly, the incidence matrix for Ri can be derived from 
the incidence matrix of X. 

The incidence matrix of A = Xi U Xi + \ appears in the upper left corner of the matrix for R := Ri. 
We obtain the matrices for the linear maps aj : H r (X ti ) — > H r (X Si ) and /3j : H r (X ti+1 ) — > H r (X Si ) by 
using the persistence algorithm Q [19] on R and A as follows. First, we run the persistence algorithm on 
the incidence matrix for A to compute a base of the homology group H r (A). We continue the procedure 
by adding the columns and rows of the matrix for R to obtain a base of H r (R). It is straightforward to 
compute a matrix representation of the inclusion induced linear map H r (A) — > H r (R) with respect to the 
bases computed by the persistence algorithm. 

Step 2. Step 2 takes the matrix representation M Pr constructed out of matrices «j , /3j computed in step 
1, and uses four elementary transformations Ti(i),T2(i), T^(i), and T^i) defined below to transform M Pr 
to M p i = T...(- • • )M Pr , whose total number of rows and columns is strictly smaller than that of M Pr . For 
convenience, let us write p = p r and p' = p' r . Each elementary transformation T modifies the representation 
p to the representation p' while keeping indecomposable Jordan cell representations unaffected but possibly 
changing the bar code representations. Some of these bar code representations remain the same, some are 
eliminated, and some are shortened by one unit as described below. For each elementary transformation we 
record the changes to reconstruct the original bar codes. The elementary transformations are applied as long 
as the linear maps or satisfy some injectivity and surjectivity property. When no such transformation 
is applicable, the algorithm terminates with all a,i and fa being necessarily invertible matrices. At this point 
the bar codes can be reconstructed reading backwards the eliminations/modifications performed. The Jordan 
cells then can be obtained as detailed in Step 3. 

The elementary transformations modify the bar codes as follows: 

• Ti(i) shortens the bar codes(i — 1, k} to (i, k} if i > 2 and shortens the bar codes (m, k}, m < k, to 
(1, k — m} if i = 1. 

• T2{i) shortens the bar codes {I, i + km] to {I, i — 1 + km] for k > 0. 

• Tz(i) shortens the bar codes [i, k} to [i + 1, k} for i < m and to [1, k — m} if i = m. 

• Ti(i) shortens the bar codes {I, (i + 1) + km) to {I, i + km) for k > 0. 
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It is understood that if an elementary transformation applied to a bar code provides an interval which is not a 
bar code, then the bar code is eliminated. Consequently T\(i) eliminates the bar codes (i — 1, i), (i — 1, ij^, 
T2(i) eliminates the bar codes (i — 1, £], T^(i) eliminates the bar codes + 1), and T±(i) 
eliminates the bar codes (i, i + 1), [i, i + 1). 

To decide how many bar codes are eliminated one uses Proposition 16. 1 1 below. Let #{i,j} p denote the 
number of bar codes of type {i, j}. 



Proposition 6.1 

1. i + l) p = dimker /3j H ker aj + i 

2. #M P = dim(V 2i /m(V 2i+1 ) + aiQfa-i)) 

3. + 1} P = dim(^(V 2i+1 ) + ^(kerft-i)) - dim^ i (V 2i+1 )) 

4. i + l) p = dim(a i (V r 2 j-i) + ft(ker a i+ i)) - dim(a»(V2t-i)) 

The following diagrams define the elementary transformations and indicate the relation between the rep- 
resentation p = {Vi, ctj, Pi} and the representation p' = {V(, a[, (}[} obtained after applying an elementary 
transformation. 



• Transformation T\(i): 



"i+i T/ ft at T7 ft-i 
6 ^21+ 1 V2i ■< V-?, 



'2i-l 



V 2i - 



Vii-i = V 2i -i/ ker(A-i), ^ = ^/^(ker^-i)), V fc = V fc ', A: ^ 2i, 2i - 1 



Transformation T 2 {i): 

■ ■ ■ ^ v 2i+1 v 2i ^a- v 2l . x v 2l _ 2 



K 2i 1/ 2i-l 



V 2 'i = Pi{V 2i +i), V 2 \_ x = ar^ft^M+i)), Vfc = V£, fc ^ 2i - 1, 2i 



Transformation Ta(i): 

ft+l t r a i+l T/ ft t r Oii ft-1 

• • • *" V2i+2 "« V2i+1 V 2 i ■< V 2 i-l > 



ft 

v 2i+l ^ ,/ 2i 



V^'j = aiWn-i), V> i+1 = pr^OiiV^i)), V k = V£, k^2i,2i + l 



if i — 1 eliminates the bar codes (m, m + 1) and (m, m + 1] 
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Transformation T±{i): 

Pi+l - ff a i+l T/ ft T , Qj T , 

• • • V2i+2 ■* V2i+1 V2i ~*- V2i-1 s 




Pi 



VL+i = V 2i+i/ = F 2i /ft(ker(a i+ i)), V k = V>, k / 2i, 2i + 1. 



The verification of the properties stated above and the proof of Proposition 16.11 are straightforward for 
indecomposable representations described in section |4]and therefore for arbitrary representations. 

As one can see from the diagrams above, when ft_i is injective, the representations p and p' are the 
same and we say that T\{%) is not applicable. Similarly, when ft is suijective, T 2 (i) is not applicable, when 
cti is surjective, T^{i) is not applicable, and when ctj+i is injective, T^ii) is not applicable. When all aj, ft 
are invertible, no elementary transformation is applicable and at this stage the algorithm (Step 2) terminates. 

To explain how the algorithm works, it is convenient to consider the following block matrices B 2 j_i and 



B 2 i, i = 1, 
-ft. Let 



, m 



, which become the sub-matrices of M Pr in ( [3]) when the entries ft are replaced with 



D 



2i-l 



for i = 1, 2, • • • (to — 1) and 



2m- 1 



«i ft 

Oi+i 



Q^m ftn 

«i 



B 



2i 



B. 



2m 



ft o 

ft+i 



/3m 

ai ft 



(V) 



(8) 



We modify M p by modifying successively each block B^. When m > 1 the algorithm iterates over the 
blocks in multiple passes. In a single pass, it processes the blocks Bi, B 2 , . . . , B 2m in this order. 

WhenB^.D- 



ft;, 



ft 



is processed then: 



1. If ft_i is not injective, we apply T\(i). This boils down to changing the bases of Vn-\ and Vn so that 
the matrix B 2 (j_i) becomes 



#-1,1 










«i,2 


Pi 







ft. 



with (ft-i,i 0) in column echelon form and 



(.2 





in row echelon form. 



In this block matrix the first and third columns correspond to V^ i _ l and V^+i respectively, and the 
first and third rows to V2(i_i) and respectively. The second column and row become "irrelevant" 

'P'i-i Q\_fPi-i,i o 



as a result of which the modified block matrix B 2 (j_i) becomes 



P'i 



Pi 



2. If ft is not surjective, we apply T 2 (z). This boils down to changing the bases of Vn-\ and V^i so that 



the matrix B 2 ( i _ 1 ) becomes 



f_Pi-i,i 

V 



ft 



ft 



8,1 







i-1,2 



ft 



i,2 
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with ~\ in row echelon form and (a? x 0) in column echelon form. 

In this block matrix the second and third columns correspond to V^i-x and V^+i respectively, and the 
first and second rows to V2U-1) and respectively. We make the first column and third row "irrelevant" 



as a result of which the modified block matrix £?2(i-i) 
When B21-1 is processed then: 



becomes 



O ; 



«i,2 



3. If a, is not surjective, we apply T^{i). This boils down to changing the bases of V21+1 and V21 so that 
the matrix B2%-\ becomes 





a 1 


#1 


#,2 















tti+1,1 





with KqM m row echelon form and (Pf^ 0) in column echelon form. 



In this block matrix the first and third columns correspond to V-a-i and V^ i+1 respectively, and the first 
and third rows to and V^+2 respectively. We make the second column and second row "irrelevant" 



as a result of which the modified block matrix i?2i— 1 becomes 



d Pi \ (<*\ # 2 







a i+i 



aj+1,2 



4. If is not injective, we apply T^i). This boils down to changing the bases of V^+i and V21 so 
that the matrix .821-1 becomes 





0i,i 


Ph 




Pi 








oti+1,1 






with 0) in column echelon form and 



Ph 




in row echelon form. 



In this block matrix first and second columns correspond to V21-1 and V^ i+1 respectively, and second 
and third rows to and V 2 (j + i) respectively. We make the third column and first row "irrelevant" as 



a result of which the modified block matrix B 



2i-l 



becomes 



«J Pi\ _ H Pi 



a 



i+l 







Explicit formulae for a's and [3's are given at the end of this section. At each pass the algorithm may 
eliminate or change bar codes, and if this happens, the matrix has less columns or rows. If this does not 
happen, the algorithm terminates, and indicates that there is no more bar code left. At termination, all a, 
and become isomorphisms. The bar codes can be recovered by keeping track of all eliminations of the 
bar codes after each elementary transformation. A bar code which is not eliminated in a pass gets shrunk by 
exactly two units, during that pass, that is, a bar code {i,j} shrinks to {i + 1, j — 1} by exactly two distinct 
elementary transformations, by elementary transformations. For example if m = 5 the bar code (1,5] 
during the pass became (2,4] as result of applying T\(l) when inspecting B\ and T 2 (5) when inspecting 
B 9 . 

When a bar code [i, i\ is eliminated, say, in the /cth pass, we know that it corresponds to a bar code 
[i — k + 1, i + k — 1] in the original representation. Similarly, other bar codes of type {i, i + l} eliminated 
at the kth pass correspond to the bar code {i — k + l,i + k}. In both cases, the multiplicity of the bar codes 
can be determined from the multiplicity of the eliminated bar codes thanks to Proposition 16.11 
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When m = 1, the operations on above minors are not well defined. In this case we extend the quiver 
G2 to Ga (m = 2) by adding fake levels *2,S2 where H r (X t2 ) = H r (X S2 ) = H r (X Sl ) and 012, P2 are 
identities! 



A high level pseudocode for the step 2 can be written as follows: 

Algorithm BarCode(M p ) 

Consider the block sub-matrices B\, . . . , B m of M p ; 
Repeat 

for j := 1 to 2m do 

1. if j = 2i — 1 is odd 

A. if oij+i is not injective, update 7?2i-i := Td(i)(B2i-i). 

B. if «j is not surjective, update i?2i-i := T^{i){B2i-\). 

C. delete any rows and columns rendered irrelevant. 

2. if j = 2i is even 

A. if /3j + i is not surjective, update B2i := T2(i)(B2i). 

B. if /3j is not injective, update B21 '■= T\(i)(B2i)- 

C. delete any rows and columns rendered irrelevant, 
endfor 

until M p is not empty or has not been updated. 
Output M p . 

Example. To illustrate how step 2 works, we consider a representation given by 



The reader can notice that this is the representation p\ for a simplified version of the example provided in 
Fig 2 with the cylinder between the critical values 62 and #3 removed. 

• Inspect Bi and B2 . No changes are necessary. 

• Inspect -B3. Since as is not injective , one modifies the block by applying T^{2) which makes both 



• Inspect the blocks Ba, B§, Bq, Bj. No changes are necessary. 

• Inspect B%. Since f5± is not injective, one modifies the block by applying 71(1) which leads to a\ = 




(9) 



and fa equal to 





7 Other easier methods can also be used in this case 
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Indeed the block B$ is given by 



Ik 







ai I3i 



( 


1 











^ 







1 













1 


1 


2 


1 







-3 


4 


2 





1 


V 


-2 


1 


2 





o J 



Since p± is already in column echelon form one only has to change the base of V 2 to bring the last 
column of a 1 in row echelon form which ends up with 



Therefore a[ 



-4 3 
-3 



1 











\ 





1 











1 


1 


2 


1 





-4 


3 





-1 


1 


-3 








-1 


/ 








<-\ 






& 


= 1 







The algorithm stops as all a[s and /^'s are at this time invertible. The last transformation ?i(l) has elim- 
inated only the bar code (4, 5], and the previous, which was the first transformation, T±(2), has eliminated 
only the bar code (2, 3). This can be concluded from Proposition 16. II In view of the properties of these two 
transformations, one concludes that these were the only two bar codes. 



Step 3. At termination, all m and become isomorphisms because otherwise one of the transformations 
would be applicable. The Jordan cells can be recovered from the Jordan decomposition of the matrix 

-1 „ a-l a-l „ a-l „ o-l „ a-l 



T = & , ■ Q!j_l • Pi 



Pi 1 ■ 011 ■ Pn 



Oi r 



■p 



1 

i+1 



o 



1+1 



■Pi • ol\ for any i. 



Standard linear algebra routines permit the calculation of the Jordan cells for familiar algebraic closed fields. 
Note that if k is not algebraically closed, Step 1 and Step 2 can still be performed and the matrix T can be 
obtained. In this case it may not be possible to decompose the matrix T in Jordan cells unless we consider 
the algebraic closure of k. It is however possible to decompose the matrix T up to conjugacy as a sum of 
indecomposable invertible matrices while remaining in the class of matrices with coefficients in the field k. 
This is the case for the field k = Z2. 



In the Example above T 



provides the Jordan cell (A = 3,k = 2) 



6.1 Implementation of Ti(i),T 2 (i), T 3 (i) and T 4 (i). 

'Pi-i 



1. T\(i) acts on the block matrix -82(1-1) = ( ^' 1 X ) . First we modify -82(1-1) to the block matrix 

\ a i Pi J 

° ! J where 0) = Pi-i ■ R(Pi-i) and (a^i a^ 2 ) = «i • R(Pi-i). Recall 

the definition of R(-) and L(-) given under notations in the introduction. Then, one passes to the block 
matrix 



a 



1,1 







a 



i.2 

2 

i.2 



a 



i.2 



Pj with 



PI 



( "l 1 ^ = L(a ij2 )-a it i and f^ 2 ) 



a 



i.2 
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The modified block matrix is y 1 o 1 ' 1 | . 

V °£i P* J 

2. T2(i) acts on the block matrix B 2 u-i\ = y 1 " 1 ® J . First we modify £?2(i-i) to the block matrix 



^ V # j where = L(&) ■ A and Q 

Then, one passes to the block matrix 



' P' 11 P' 12 \ 

4 1 ' # J with ("j 1 ) = a^^Ki), Q^) = Oi.ii?^!), and = /^(a^) 



"!,1 «i,2 



a?! 
The modified block matrix is ( ^"j 1 ' 2 ^ 

3. Tafi) acts on the block matrix B 2 %-x = \2 ^ V First we modify B^i-x to the block matrix 

\ a i+1 J 

3? J where ^ = a ■ R{a t ) and = ft • fl(oi). 

i+l/ 

Then, one passes to the block matrix 

with {(31 0) = /3 2 ■ R((3f), {f3l /3J 2 ) = A 1 • 

and (aj+1,1 0^+1,2) = ctj+i • R{Pf)- The modified block matrix is ^ l)2 J . 

4. T4(i) acts on the block matrix Bn-x = \ 2 % ) ■ First one modifies Bn-\ to the block matrix 

Co af+i\ ^o' 2 ) wnere °) = ■ ^(«i+i) and (ft.i ft i2 ) = A ■ R(a i+ x). 

Then, one passes to the block matrix 

2 P J with (^ 2 ) = L(A, 2 )-ft,2, (^J = i(A,2)-/3i,i and = L(A |2 )-«i- 

a i+l,l / 



# 
ai + xj 



Plx Pi 


Oi+1,1 a.j+1,2 



<X 2 



. «i+l,l 



/a 2 /? 2 N 
The modified block matrix is * h 

V «i+l,l, 
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6.2 Time complexity 

Let the input complex X have n simplices in total on which the circle- valued map / is defined which has m 
critical values. 

Then, step 1 takes 0{nd) time to detect all the critical values where d < n is the maximum degree of 
any vertex. The critical values can be computed by looking at the simplices adjacent to each of the vertices. 
To compute the matrices a,i and we set up the matrices of size 0(n) x 0(n) and run persistence on them. 
Using the algorithm of [161, this can be achieved in 0(M(n)) time where M(n) is the time complexity of 
multiplying two n x n matrice^l. Since we perform this operations for each of the critical levels and the 
spaces between them, we have 0(mM(n)) total time complexity for step 1. 

In step 2, we process the matrix M Pr iteratively until all BarCode representations are removed. In each 
pass except the last one, we are guaranteed to shrink a bar code by at least one unit. Therefore, the total 
number of passes is bounded from above by the total length of all bar codes. Theorem 13. 1 1 implies that a 
bar code cannot come back to the same level more than max Sj dim H r (X Si ) times which can be at most 
0(n). Therefore, any bar code has a length of at most 0(nm) giving a total length of 0(n 2 m) over all bar 
codes. Hence, the repeat loop in the algorithm BarCode cannot have more that 0(n 2 m) iterations. In each 
iteration, we reduce the block matrices each of which can be done with 0(M (n)) matrix multiplication time 
|[T5l . Since there are at most 0(m) block matrices to be considered, we have 0(mM(n)) time per iteration 
giving a total of 0(n 2 m 2 M(n)) time for step 2. 

Step 3 is performed on the resulting matrix from step 2 which has 0(mn) x 0(mn) size. This can again 
be performed by matrix multiplication which takes 0(M(mn)) time. 

Therefore, the entire algorithm has time complexity of 0(m 2 n 2 M(n) + M{mn)). 

7 Conclusions 

We have analyzed circle valued maps from the perspective of topological persistence. We show that the 
notion of persistence for such maps incorporate an invariant that is not encountered in persistence studied 
erstwhile. Our results also shed lights on computing homology vector spaces and other topological invari- 
ants from bar codes and Jordan cells (Theorems 13.1 l and I3T21 . We have given an algorithm to compute the 
bar codes and the Jordan cells; the algorithms can also be adapted to compute zigzag persistence. In a 
subsequent work, Burghelea and Haller have derived more subtle topological invariants like Novikov ho- 
mology, monodromy Q, Reidemeister torsion, and others from bar codes and Jordan cells confirming their 
mathematical relevance. We have not treated in this paper the stability of the invariants; see for partial 
answer. 

The standard persistence is related to Morse theory. In a similar vein, the persistence for circle valued 
map is related to Morse Novikov theoiy lTT8l . The work of Burghelea and Haller applies Morse Novikov 
theory to instantons and closed trajectories for vector field with Lyapunov closed one form J2]. The results 
in this paper will very likely provide additional insight on the dynamics of these vector fields and have 
implications in computational topology in particular and algebraic topology in general. 
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Appendix 

In this Appendix we explain the calculation of the r-invariants for the example depicted in Fig 2. The 
representation po has vector spaces that are all one dimensional and maps ctj = that are all identity. 
Hence, there is no bar code, but one Jordan cell A = 1, k = 1. 

It is not hard to recognize from Fig 2 that the maps for the representation p± are given by: 















[-1 


4 


V 


| ; «2 = 









1 


V 




v0 


0/ 



01 





1 

1 

1 
1 



1 

1 

1 
1 



We proceed with the step 2 of the algorithm. 

• inspect B\ - no change for p = pi; inspect B2- no change. 



inspect B3, - since 02 is not surjective apply 7/3(2). This changes ct2,/32,ot3 into a' 2 

,1 o^ 

^2 = ( n 1 ) ' a 3 = I 1 I • Update and continue. 
0, 



1 
1 



,° 1 

inspect B4- no changes. 



inspect B§ - since 03 is not surjective, apply T^(3). This changes and (3^ into a' 3 = M and 



f3' 3 = (^Q -j* J ■ Update and continue. 
• inspect Bq- no changes. 
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• inspect B-? - since «5 is not injective, apply T±(&). This changes /?4 and 05 into c/ 5 



(■ 



1 
1 



') 



and 



1 
1 



. Update and continue. 



• inspect Bg - no change; inspect Bg- no change; inspect B\q- no change; inspect B\\- no change. 



Since at this time all a-s and /3-s are invertible, step 2 terminates. 

Book keeping. The last transformation Ti(l) has eliminated the bar code (9q, 6\ + 2ir] (by Proposition 
16.11 ) and nothing else. This bar code was not the modification of any other bar code by the previous elemen- 
tary transformations. The previous transformation T±(&) has eliminated the bar code (64,65) and nothing 
else (by Proposition 16. II ). This bar code was not the modification of any other bar code by the previous 
transformations. The transformation T3(3)has eliminated the bar code [#3,6*3] (by Proposition 16.11 ) which 
was the modification of [#2 ,#3] by T^(2). These are all bar codes as listed in the table in section [3] To 

calculate the Jordan cells we use step 3. We calculate the Jordan cells of f ^ jjj • ( f j J J ) _1 which 

is (A = 3, k = 2) as listed in the table in section|3] 



• inspect B12 - since (3q is not injective, apply This changes /3q, a±, fi\ to f3' e 
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